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abstract 



In this work we present a mathematical framework that quantifies the gas- 
exchange processes in the lung. The theory is based on the solution of the one- 
dimensional diffusion equation on a simplified model of lung septum. Gases 
dissolved into different compartments of the lung are all treated separately 
with physiologically important parameters. The model can be applied in mag- 
netic resonance of hyperpolarized xenon for quantification of lung parameters 
such as surface-to-volume ratio and the air-blood barrier thickness. In gen- 
eral this model provides a description of a broad range of biological exchange 
processes that are driven by diffusion. 
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Introduction 



Gas exchange is the essential process that happens in the lung. However, 
quantifying gas exchange has not been an easy task, mainly because of the 
complicated structure of the lung. In general, the human lung can be seen 
as a porous media with capillary blood flows inside, which mainly consists of 
the red blood cells (RBCs) and plasma. The blood flows are separated from 
the air space by several layers of tissue, including epithelium, endothelium 
and interstitium [1]. These tissues together are called the air-blood barrier. At 
equilibrium gas molecules in the alveolar space are under constant exchange 
with those dissolved into the air-blood barrier and blood in the lung. 

Gas-exchange is the result of diffusion. For most gas species their concentra- 
tions are much higher in the alveoli than in the parenchyma they dissolve 
into (depending on the corresponding Ostwald solubility), and equilibrium is 
established once the tissue and blood are saturated with the gases. In many sit- 
uations, due to the different chemical environment the dissolved gas molecules 
experience, they can usually be distinguished from the free gas, which allows 
selective manipulation (e.g., RF saturation) of either state of the gas. Such 
manipulations usually lead to a broken equilibrium between the dissolved and 
free gases, and because gas exchange will eventually re-establish the equilib- 
rium, one can use this technique to quantify the gas-exchange processes [2]. 

The current work is largely motivated by the recent development of mag- 
netic resonance (MR) of hyperpolarized 129 Xe(HXe) [3]. HXe makes an ideal 
contrast agent for quantifying gas exchange in the lung because of its large 
chemical shift when dissolved into lung tissue and blood plasma, both at 197 
ppm[4] relative to the frequency of free xenon. In blood, xenon also dissolves 
into the RBCs and binds hemoglobin, which gives rise to yet another unique 
peak at 217 ppm [4]. A simplified model of lung septum with the source of 
each dissolved-xenon peak identified is shown in Fig. 1. Because of the large 
chemical shifts of the dissolved xenon in the lung, one can use the technique of 
chemical shift saturation recovery (CSSR) [5] to selectively saturate the dis- 
solved xenon and then monitor the recovery of the xenon signals at both 197 
and 217 ppm as functions of exchange time, which we call the xenon uptake 
dynamics. 

The growth curves of the dissolved-xenon signals carry important information 
of the lung. Several groups have developed theories of the uptake dynamics 
using important lung parameters [6,7]. However, due to either limited range of 
exchange times or limited field strengths that were unable to distinguish the 
two peaks of dissolved xenon, the previously developed theories were not fully 
tested against the uptake dynamics of both dissolved-xenon compartments. 
In this work we presented a general theory of gas exchange in the lung with 
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a greatly simplified geometry of lung septum. This theory first calculates the 
dissolved xenon in the air-blood barrier and in the blood separately; then, in 
order for it to be used for xenon uptake dynamics, the xenon in blood is further 
separated into the RBC xenon and plasma xenon. Finally the barrier (tissue) 
xenon and the plasma xenon are combined for their common chemical shifts. 
We note that although this theory is aimed for interpreting gas exchange in 
the lung, the method can also be applied to general exchange processes in a 
biological system driven by diffusion. 



Theory 



We begin with the 1-D diffusion equation within a simplified model of lung 
septum of thickness d, as shown in Fig. 1. If we consider all the gas dissolved 
into lung tissue as a single component, i.e., we don't distinguish between gas 
in tissue and blood, then the dissolved gas density, M d , as a function of the 
distance x from the edge of the tissue and time t, satisfies the diffusion equation 

dM d=D d 2 M d 



dt dx 2 ' 

subject to the initial condition 

M d (x,0) = ,xe(0,d), (2) 

and the boundary conditions 

M d (0, t) — M d (d, t) — AMf . (3) 

In the above equations, Mf is the density of free gas (at ppm), D is the 
diffusion coefficient of dissolved gas, d is the septal thickness, and A is the 
Ostwald solubility of xenon in lung septum. We also made the assumption 
that all the dissolved gas molecules that hit the wall at a certain time t leave 
the tissue (i.e., exchange with the free gas molecules). 

The complete solution of this problem is well known and can be represented 
by a sum of infinite series: 

x ,„f 4 ^ 1 (2n-l)nx p ^-if^ A 
M d (x, t) = AM f 1 - - V sin ^ —!— e' D —^ 1 . (4) 



d 



In a real experiment the density distribution M d (x,t) can be converted in the 
signal distribution S d (x,t) of the dissolved gas. S d (x,t) is proportional to the 
total surface area Sa in the lung, and can be normalized by the total free-gas 
density, MfV g , where V g is the gas volume in the lung. If we define the gas 



4 



tissue 

(197 ppm) 






Red blood cells 

(217 ppm) 

and 
blood plasma 

(197 ppm) 





d-d 



d 2 



Fig. 1. A simplified model of lung septum. 
exchange time constant T as 

T = 

then Eq. (4) can be simplified to 
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We now consider the dissolved-gas signal from each compartment in the lung. 
As shown in Fig. 1, we use a simple model that assumes a layer of tissue 
(air-blood barrier) of thickness 8 (8 < d/2) at each side of the septum. The 
total signal from the tissue, denoted by Sdi(t), can be calculated as the spatial 
integral of Sd in Eq. (6) over the two regions from to 5 and from d — 8 to d. 
If we let fi = XdS^/Vg and k — 8/d, then the the gas signal from the tissue 
compartment S^i is 



Sdi(t) = v{ 2« 



E 
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(7) 



The gas signal from the blood is more difficult to calculate due to the flow 
effect. Assuming no flow, the gas signal from the static blood, Sd2s(t), is simply 
the spatial integral of Sd over the region from 5 to d — 5: 



Sd2s(t) 
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To deal with the partial volume effect caused by the blood flow, we follow the 
method given by Patz in [7]. In brief, we divide the region with dissolved gas 
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Fig. 2. Lung septum with blood flow. 

into three sections, as shown in Fig. 2. First we introduce a new quantity, the 
pulmonary capillary transit time, defined as the average time an RBC spends 
in the gas-exchange zone (i.e., in contact with the alveolar space shown in Fig. 
2) in the lung, and denoted by tx- If we let t be the exchange time, and r be the 
time a certain infinitesimally thin layer of blood spends in the gas-exchange 
zone, then the signal contribution from the blood gas in this thin layer is, 
in terms of Sd2s, S^irfdr/tx, and therefore the total contribution of this 
"partial volume" (the two green areas in Fig. 2) can be calculated with a time 
integral over Sd2s; on ^ ne other hand, for t < tx, the region of (1 — t/tx) thick 
(blue area in Fig. 2) is always in contact with the gas-exchange zone during 
t, thus the signal with this section can still be treated as static. Therefore the 
total signal is 



S d2 = 2 f ^S d2s (r) + S d2s (t) 
■to *x \ t x J 

E —± cos(rwrK) (l 
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where I have used the integral 



' e -n^/T dT ^T e _ nH/T 

o n 2 



(10) 



In order for this theory to be applied to dissolved-HXe uptake dynamics we 
need to find expressions for the signal amplitudes at 197 ppm and 217 ppm. 
As shown in Fig. 1, for humans, the HXe in tissue and blood plasma share the 
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same resonant frequency at 197 ppm, and HXe in the RBCs resonates at 217 
ppm. Thus if we let r\ be the fraction of dissolved gas in the red blood cells, 
the peak amplitude at 197 ppm, S^p, combining HXe in tissue and plasma, is 



S T p(t) = S dl (t) + (l-r])S d2 (t) 



= A^-^- 2 E f . 

k n=odd 




+ 2//(l 



(11) 



and the peak amplitude of HXe in the RBCs at 217 ppm, Srbc, is 



Srbc(*) = VSd2(t) 




(12) 



Discussion 

Test of the theory using HXe experiments will be presented separately. Nev- 
ertheless, insights can be gained by plotting 5rp and Srbc i n Eqs. (H) an d 
(12), respectively, as functions of the gas-exchange time for n < 9, as shown 
in Fig. 3. The plots are made using the following estimated parameters for the 
human lung: fi = 0.03, k = 0.15, rj = 0.5, T = 30 ms and tx — 1-5 s. These 
two curves, by comparison, are similar in shapes and characteristics with the 
previously published dissolved-HXe data [6,8,9], which supports the validity 
of the presented theory. We would like to pointed out two interesting features 
of the curves. First, for exchange time longer than 100 ms, both curves grow 
linearly with time. This is because both tissue and blood are saturated with 
gas (HXe in this case) at about 100 ms and the signal growth after which is 
purely due to blood flow; second, although blood is never is directly contact 
with the free gas in alveoli, there is still signal from RBC at the very short of 
exchange time. This is a result of diffusion — the propagation of probability 
density in diffusion is non-local. 

The presented theory can potentially be used to fit the experimental data 
and extract important parameters related to lung physiology and function. 
The parameters /z, k and tx represent, respectively, the product of septal 
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Fig. 3. Plots of Stp and 5rbc hi Eqs. (11) and (12), respectively, for n < 9, as 
functions of the gas-exchange time up to 700 ms. The parameters used are: fj, = 0.03, 
k = 0.15, rj = 0.5, T = 30 ms and tx = 1-5 s. The insert show details of the curves 
at short exchange times. 

thickness d and the surface-to-volume ratio Sa/V^, the ratio between the air- 
blood barrier 5 and septal d thicknesses, and the pulmonary capillary transit 
time. 7] is the partition of the RBC HXe in the blood, which can be used to 
calculate the hematocrit (Hct) using 

Hct = — — — — , (13) 

WArbc + (1 - V)/ x p 

where Arbc an d Ap are the Ostwald solubilities of HXe in the RBCs and 
plasma, respectively. 

The current model used a fair number of simplifications, including the naive 
geometry of the septum, undistinguished gas solubilities and diffusion coeffi- 
cients in the tissue and blood. Therefore, great improvements can be done by 
carefully dealing with these factors. 



Conclusions 



We have developed a theory of gas-exchange in the lung based on a simple 
model of 1-D gas diffusion. This theory carefully treats the dissolved gas in 
different compartments of lung septum with parameters of physiological im- 
portance. It can be used to quantify pulmonary function using the dynamics 
of the dissolved hyperpolarized xenon in the lung. 
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